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Recent studies show that the D" layer, just above the Earth's core-mantle boundary, is 
composed of MgSi0 3 post-perovskite and has significant lateral inhomogeneity. Here we 
consider the D" diversity as related to the single-crystal elasticity of the post-perovskite 
phase. We measure the single-crystal elasticity of the perovskite Pbnm-CaW0 3 and post- 
perovskite Cmcm-Calr0 3 using inelastic X-ray scattering. These materials are structural 
analogues to same phases of MgSi0 3 . Our results show that Cmcm-Calr0 3 is much more 
elastically anisotropic than Pbnm-Calr0 3 , which offers an explanation for the enigmatic 
seismic wave velocity jump at the D" discontinuity. Considering the relation between lattice 
preferred orientation and seismic anisotropy in the D" layer, we suggest that the c axis of 
post-perovskite MgSi0 3 aligns vertically beneath the Circum-Pacific rim, and the b axis 
vertically beneath the Central Pacific. 
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Since the discovery of the transition from perovskite (Pv) to 
post-perovskite (pPv) in MgSi0 3 (refs 1,2), pPv-MgSi0 3 has 
been considered to be a major component of the D" layer 
just above the core-mantle boundary (CMB). Seismological 
observations of the D" layer are difficult to interpret as they 
suggest both non-uniform response at the layer boundary and 
non-uniform anisotropy inside the layer 3-5 . With this background, 
elastic anisotropy or single-crystal elasticity of pPv-MgSi0 3 
has been the focus of theoretical calculations 6-10 . However, no 
experimental data are available on the single- crystal elasticity of 
pPv-MgSi0 3 , because it is unstable at ambient pressure. 

Cracra-CaIr0 3 or pPv-CaIr0 3 has been frequently selected as a 
representative analogue of pPv-MgSi0 3 (refs 11-13). According 
to the phase diagram 14 , pPv-CaIr0 3 is the stable phase at ambient 
conditions and Pv-Calr0 3 is the high- temperature and low- 
pressure phase. Here we report the first experimental data on 
single- crystal elasticity of Pv and pPv structures in a 
homogeneous chemical composition of Calr0 3 . 



Results 

Characterization of single-crystal specimen. Figure la,b shows 
single crystals of Pv-Calr0 3 and pPv-CaIr0 3 used in the present 
study. The Pv-Calr0 3 sample was synthesized at 2GPa and 
1,450 °C for 15 h in a piston cylinder apparatus at Okayama 
University. We produced many grains with size of a few tens of 
microns, although most of them turned out to have multiple 
domains. We examined more than 20 grains using a four- circle 
X-ray diffractometer at Yamaguchi University and selected a 
quasi-single-domain crystal with size of ~ 20 urn. Its lattice 
constants were determined as a = 5.3527(1), fr = 5.5969(5) and 
c = 7.6804(6) A, which yields density (p) of 8,091 kgm -3 . The 
pPv-CaIr0 3 sample was synthesized by slow cooling from 
1,000 °C in CaCl 2 flux 15 . It had a pine-needle shape with length of 



a few 100 microns and width of a few 10s of microns. Its lattice 
constant and densities were respectively a = 3.145(2), 
b = 9.861(6), and c = 7.297(5) A, and p = 8,211 kg m" 3 , which 
are consistent with literature values 16 . 



Inelastic X-ray scattering measurement and results. Figure lc,d 
compare experimental velocities obtained from the inelastic X-ray 
scattering (IXS) spectra with velocities calculated from the opti- 
mized elastic constants. The details of the IXS measurement are 
described in the Methods section for the measurement of IXS at 
SPring-8 BL35XU; one of the important advantages of IXS 
compared with Brillouin scattering is its geometrical freedom 
(see Supplementary Methods). The analytical procedures for 
determining single- crystal elastic constants are described in 
the Methods section for data analysis for elastic constants. 
The measurement data are summarized in Supplementary 
Tables 1 and 2. 

Table 1 summarizes the results of elastic constants and gives 
isotropic averages. The values are reasonably consistent with the 
theoretical results 17 except for a few discrepancies. The high- 
pressure phase of pPv-CaIr0 3 is less compressible than Pv- 
Calr0 3 whereas the theoretical prediction suggests the opposite 
relation. In addition, isotropic averages of V P and V s are 
significantly larger for pPv-CaIr0 3 than for Pv-Calr0 3 
according to the present results. Figure 2 compares the present 
results and theoretical calculations for the velocity surfaces of Pv- 
Calr0 3 and pPv-CaIr0 3 . We see fairly good agreement between 
the two results for pPv-CaIr0 3 , although we recognize significant 
differences between them for Pv-Calr0 3 . Figure 3 shows the 
compatibility of the elastic constants with compression data 18 . 
The present data are consistent with the compression data for Pv- 
and pPv-CaIr0 3 ; the theoretical result for Pv-Calr0 3 , however, 
yields a different slope from the others. 
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Figure 1 | Specimen and measurement data. Left and right columns refer to Pbnm- and Cmcm-Calr0 3 , respectively. (a,b) Photographs of single crystals 
used. The horizontal edge lengths are ~5mm and ~800|im, respectively, for a and b. (c,d) Overall picture for the experimental data and analytical 
fitting. Black lines indicate the direction of wave propagation in the crystal, and blue triangles and red circles are experimental velocities and recalculated 
values after least-squares analysis. The vertical direction is positive z direction, while horizontal leftward and rightward directions are positive x and y 
directions, respectively. Note that the difference between experimental and recalculated velocities is less than 5%. The misfit uniformly distributed 
indicates untwinned crystal specimens. 



2 



NATURE COMMUNICATIONS | 5:3453 | DOI: 10.1038/ncomms4453 | www.nature.com/naturecommunications 
© 2014 Macmillan Publishers Limited. All rights reserved. 



NATURE COMMUNICATIONS | DPI: 10.1038/ncomms4453 



ARTICLE 



Table 1 | Single-crystal elastic moduli of Pbnm- and Cmcm-Ca\r0 3 . 
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Several isotropic moduli and velocities are shown on the right-hand side. The upper two rows are for Pbnm or Pv structure, while the lower two rows are for Cmcm or pPv structure. In each case, the first 
row gives the present experimental results with accuracy and the second row the theoretical results obtained by Tsuchiya and Tsuchiya 17 . The unit is 'GPa' except for velocities which have a unit of 
'ms~ 1 ' in the last two columns. The numbers in parentheses correspond to fitting error in the least square analysis. 
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Figure 2 | Velocity surface of Pbnm- and Cmcm-Calr0 3 . Upper (a) and 
lower (b) plots correspond to Pbnm- and Cmcm-Calr0 3 , respectively. Black 
lines are drawn for the present results, while blue lines are drawn for 
the results obtained by Tsuchiya and Tsuchiya 17 . Upper lines correspond 
to the P-wave velocity, and lower lines to the S-wave velocity. Open circles 
'o' and crosses ' + ' for V s specify in-plane and out-of-plane particle 
motions. 



Discussion 

Figure 4 presents stereo projections of elastic wave velocities for 
Pv and pPv structures of Calr0 3 . It shows that the ranges of 
variation of V P , the average of V s values, and the difference in V s 
values are -600, -300 and -500ms -1 for Pv-Calr0 3 and 
- 1,000, -800 and - 1,100ms" x , for pPv-CaIr0 3 , respectively. 
We see that pPv-CaIr0 3 has approximately double the anisotropy 
of Pv-Calr0 3 in terms of the elastic wave velocity. For example, 
anisotropy in the polarization shear wave velocity is —30% in 
pPv-CaIr0 3 while it is only — 15% in Pv-Calr0 3 . 

From Fig. 4d,e for pPv phase, the b axis is found to be the 
direction of the lowest Vp, while the highest Vs is between the a 
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Figure 3 | Comparison of compression curves. Upper (a) and lower 
(b) plots correspond to Pbnm- and Cmcm-Calr0 3 , respectively. The lattice 
constants a, b, and c are shown in blue, green and red respectively. Solid 
and dotted lines are drawn for the present data and Tsuchiya's theoretical 
data 17 , respectively, while symbols are the compression data. 



and c axes. These pronounced features of the pPv elasticity are 
understandable as alternate stacking of compliant and stiff 
layers perpendicular to the b axis; the stiff layer has the largest 
shear wave constant and lowest longitudinal wave constant 
between a and c axes. We recognize similarity between a and c 
axes in pPv-CaIr0 3 elasticity, while the a axis is slightly stiffer 
than the c axis in terms of both longitudinal and shear wave 
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Figure 4 | Velocity surface of Pv- and pPv-Calr0 3 . Stereo projections of elastic wave velocities for Pv- and pPv-Calr0 3 based on the present results. 
The crystallographic directions, a, b, and, c, are identical in each plot as shown in (a). The unit of scale bars is ms"l (a-c) V P , average of two V s values, 
and difference between two V s values, respectively, for Pv-Calr0 3 . (d-f) Same as a, b, and c but for pPv-Calr0 3 . 



constants. The elasticity of Pv-Calr0 3 resembles that of 
Pv-MgSi0 3 under lower-mantle conditions more than under 
ambient conditions (Fig. 5). Elasticity of Pv- and pPv-CaIr0 3 
under ambient conditions seems to be a good indicator of that of 
Pv- and pPv-MgSi0 3 under lowermost mantle conditions; a 
conversion factor of 2.2-2.4 is found between them. 

We interpret the D" diversity observed seismologically in terms 
of the single -crystal elasticity of Pv and pPv-CaIr0 3 . Significant 
anisotropy of the seismic wave velocity has been observed in the 
D" layer, while the lower mantle above the D" layer is nearly 
isotropic with a high degree of confidence 19 . Therefore, we 
plotted the velocity difference between anisotropic pPv and 
isotropic Pv both for Calr0 3 under ambient conditions and 
MgSi0 3 under mantle conditions (Fig. 6). From the plots 
for Vp and slower V s , we see that any velocity change over the 
D" discontinuity is possible both for P and S waves depending 
on the lattice preferred orientation (LPO) characteristics of the 
pPv phase. 

4 



The seismic wave phases used in the D" study (S, ScS, Sdiff) 
travel nearly horizontally in the D" layer 4 . Therefore, the two 
polarized S waves are classified as SH and SV, whose polarizations 
are in the horizontal and vertical directions, respectively. In 
seismological study, vertically transverse isotropy (VTI) and the 
tilted transverse isotropy have been assumed to explain S-wave 
polarization anisotropy in the D" layer 4 . 

According to seismological studies, the region of the Circum- 
Pacific rim is interpreted as having VTI with V SH > V sv , while the 
region of the Central Pacific has complicated VTI without any 
fixed relation between V SH and Vsv 4 . Polarization anisotropy is 
not recognized beneath the Atlantic Ocean, or the region is 
considered as having VTI with V SH « V sv . On the other hand, 
from the polarization pattern in the a-b plane of the faster V s plot 
(Figs 2b and 7), we can expect that V SH > V S y in the D" layer if 
the c axis of pPv aligns vertically. Similarly, if the b axis aligns 
vertically, the magnitude relation between V SH and V S v can be 
complicated. 
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Figure 5 | Stereo projections of elastic wave velocities for Pv- and pPv-MgSi0 3 . The crystallographic directions, a, b, and, c, are identical in each plot 
as shown in (a). The unit of scale bars is ms -1 , although a factor of 10 4 should be multiplied for V P . (a-c) V P , average of two V s values, and 
difference between two V s values, respectively, for Pv-MgSi0 3 under ambient conditions, (d-f) and (g-i) Same as a, b, and c but for Pv-MgSi0 3 and 
pPv-MgSi0 3 at 125 GPa and 2,500 K 10 , respectively. 



We can summarize the above conclusion in an alternative way: 
For the c axis vertical case, horizontal wave propagation 
directions lie in the a-b plane, for which the in-plane shear 
mode (SH) has consistently higher wave speeds (Figs 2b and 7). If 
instead the b axis is vertical, horizontal wave propagation 
directions lie in the a-c plane, for which neither the in-plane 
(SH) nor the out-of-plane (SV) is consistently faster for all such 
propagation directions (Figs 2b and 7). 

The dominant slip system of pPv-CaIr0 3 was determined as 
[100] (010) from shear deformation experiments 11 . Since the slip 
system of pPv-MgSi0 3 has not yet been determined, we assumed 
that it is similar to that of pPv-CaIr0 3 n . This assumption leads to 
LPO of pPv-MgSi0 3 , in which the b axis aligns vertically along 
the lateral flow in the D" layer. The LPO pattern is consistent 
with the complicated polarization S-wave anisotropy beneath the 
Central Pacific (Figs 2b and 7). 

Although the texture development during Pv-pPv transition in 
MgSi0 3 has not yet been clarified, we expect the c axis alignment 
of pPv-MgSi0 3 immediately after the phase transition from Pv- 
MgSi0 3 from analogy with the case for MgGe0 3 
(ref. 20). If the c axis aligns vertically, we expect VTI with 
^sh > ^sv, which is consistent with the seismic feature beneath 
the Circum-Pacific rim. 



Employing this conceptual scheme, we evaluated the feasibility 
of the LPO model of pPv in the D" layer. The volume ratio 
between (Mg, Fe)Si0 3 Pv and (Mg, Fe)0 ferropericlase (Fp) is 
0.67:0.33 for the lower mantle 21 . We also used this ratio for (Mg, 
Fe)Si0 3 pPv and Fp in the D" layer. This means that the 
anisotropy resulting from pPv LPO should be weakened by 
approximately two-thirds in the D" layer assuming isotropy of 
Fp. The observed polarization anisotropy of the S-wave velocity is 
~ 3% at most 4 . Therefore, pPv in the D" layer is expected ~4.5% 
anisotropy to account ~ 3% anisotropy observed in the D" layer. 
This requirement is easily satisfied by partial LPO of pPv 
(Table 2). 

Figure 8 summarizes the present idea for seismic anisotropy 
and LPO structure in the D" layer. The VTI of the vertical c axis 
and that of the vertical b axis can be interpreted as transformation 
LPO and deformation LPO, respectively. It is noted that isotropic 
pPv-MgSi0 3 is a reasonable interpretation for the seismic feature 
beneath the Atlantic Ocean; an alternative interpretation may be 
the VTI of the vertical a axis because of less lateral polarization 
anisotropy than in other cases (Table 2). The difference between 
beneath the Pacific and Atlantic may be attributed to differences 
in temperature, chemical composition and strain rate. The 
proposed LPO behaviour for pPv-MgSi0 3 obtained from the 
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Figure 6 | Velocity difference between anisotropic pPv and isotropic Pv. The left column presents the results for Calr0 3 in the present study, while the 
right column presents the theoretical calculation for MgSi0 3 at 125 GPa and 2,500 K 10 We recognize common characteristics among the plots for the two 
materials. Isotropic P-wave velocities of Pv are 5,774 and 13,900 ms _1 and S-wave velocities are 2,975 and 7,200 ms _1 for Calr0 3 and MgSi0 3 , 
respectively. The red arrows beside the scale bars highlight the origin of each scale bar, where the velocities of the two phases are the same. Typical 
polarizations of the S wave are shown by arrows in the S-wave plots. The crystallographic axes are identical for all plots as shown in plot a. In the stereo 
projections, the propagating waves in the a-b and c-a planes are shown along the quarter arc and the bottom edge, respectively. Note that the two 
polarization arrows in plot e between a and b axes suggest that the two S-wave velocities are nearly the same, as shown in Fig. 7. 
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Figure 7 | Velocity surfaces of pPv-MgSi0 3 at 125 GPa and 2,500 K. 

Blue lines are used to specify theoretical results similar to Fig. 2a,b. The 
upper line is the P-wave velocity and the lower lines are the S-wave 
velocities. Open circles 'o' and crosses ' + ' for V s specify in-plane and out- 
of-plane particle motions. 
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present mineral physics perspective is consistent with seismic 
observation. 

To summarize, we presented a new data set of Pv and pPv- 
Calr0 3 single-crystal elasticity by means of IXS at SPring-8. We 
used the present results to interpret seismic wave anisotropy in 
the D" layer, and proposed a model of pPv LPO in the D" layer 
consistent with seismic observations and mineral physics 
experiments on analogue materials. The present model may 
allow sophisticated discussion of global mantle convection, which 
is triggered by heating at the core-mantle boundary and 
modulated by lateral flow in the D" layer. 

Methods 

Measurement of IXS at SPring-8 BL35XU. The small size and black colour of the 
present samples limit the techniques that can be used to investigate the elastic 
properties. In particular, the crystal sizes are too small to employ GHz ultrasonic 
methods 22 or resonant ultrasound spectroscopy 23 , where both need crystals larger 
than 100 um. In addition, the crystal opacity makes it difficult to apply Brillouin 
scattering. Therefore, IXS with an X-ray beam diameter of ~ 70 um was employed. 
Each specimen was attached to glass fibres ( ~ 100 urn in diameter) using manicure 
resin and mounted on the four-circle goniometer installed at BL35XU of SPring-8 
(ref. 24). The incident X-ray beam energy was 17.794 and 21.747 keV for Pbnm- 
Calr0 3 and Cracra-Calr0 3 measurements, respectively. 

The IXS technique 25 ' 26 allows us to determine the elastic properties of a 
material by measuring the energies of acoustic longitudinal and shear phonon 
modes. When this information is combined with the precise reciprocal space 
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Table 2 | Characteristic velocities for crystals of post-perovskite structure. 
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Figure 8 | Schematic illustration of the LPO of pPv-MgSi0 3 in the D" layer. From a mineral physics viewpoint, the Pv slab (purple region) transforms to 
pPv (light blue region) and the pPv slab then moves along the path of mantle convection with deformation (light-green region). In contrast, from seismic 
observation, the region of the Circum-Pacific rim has VTI anisotropy of V SH > V S v, while the region of the Central Pacific has complicated VTI with any 
magnitude relation; that is, V SH > V sv , V S h < V sv or V SH ~ V sv . The region of the Atlantic Ocean (white circle) has VTI with V S h ~ V sv , suggesting isotropy or 
negligible LPO of pPv-MgSi0 3 . The speculated aligned axes of pPv LPO are shown by arrows for each region; the vertical c axis beneath the Circum-Pacific 
rim and the vertical b axis beneath the Central Pacific. 



position (momentum transfer) for each measurement, we can determine the sound 
velocity 27 ' 28 . If measurements are made in many different directions, the velocity 
information can be combined with the density, allowing all the elastic constants of 
a material to be determined. In the present work, the two-dimensional, 4x3, 
analyser array of BL35XU is especially advantageous as it allows access to 12 
different directions of momentum transfer simultaneously. 

The IXS is applicable on opaque materials; it is the most important advantage of 
IXS against Brillouin scattering. Further characteristics and procedures of the IXS 
measurements are summarized in Supplementary Figs 1-5. 

Data analysis for elastic constants. A least-squares fit to Christoffel equation 
was then used to determine the elastic constants 29 . Pv-Calr0 3 and pPv-CaIr0 3 
have orthorhombic symmetry and thus nine independent elastic constants. The six 



diagonal elastic constants (C u , C 22 , C 33 , C 44 , C 55 , C 66 ) is related to nearly 
independent phonon modes propagating along crystallographic axes, while the 
three off-diagonal constants (C 12 , C 23 , C 31 ) can be determined from combination of 
a few measurements. We carefully planned the measurement to acquire sufficient 
data to constrain all nine elastic constants. Using the analyser array, we obtained 
~40 and ~60 mode energies for Pv-Calr0 3 and pPv-CaIr0 3 , respectively. 
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